Exact bosonization for an interacting Fermi gas in arbitrary dimensions 
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We present an exact mapping of models of interacting fermions onto boson models. The bosons 
correspond to collective excitations in the initial fermionic models. This bosonization is applicable 
in any dimension and for any interaction between fermions. Introducing superfields we derive a field 
theory that may serve as a new way of analytical study. We show schematically how the mapping 
can be used for Monte Carlo calculations and argue that it should be free from the sign problem. 
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The study of interacting fermionic systems in cases 
when the Landau Fermi liquid theory fails to describe 
all interesting effects is an open problem of condensed 
matter theory. Very often conventional methods [l[ are 
not efficient due to divergencies in perturbation expan- 
sions leading to the re-summation of complicated series. 

It is difficult to list here all the problems encountered 
in the study of e.g., strongly correlated systems. The 
most clear examples are provided by one dimensional 
(ID) systems where perturbative methods are especially 
inconvenient but models suggested for describing high Tc 
superconductivity, see, e.g. [2], are not simpler. Generi- 
cally, the low temperature physics of systems of interact- 
ing fermions is naturally described in terms of bosonic 
collective excitations that can be expressed only by infi- 
nite series of conventional diagrams. 

The numerical study of fermionic systems encounters 
difficulties as well. The powerful Monte Carlo (MC) 
method suffers the well known sign problem 0, 0, H, @, 01 
leading to a drastic increase in the computing time. 

All these examples call for a reformulation of inter- 
acting fermions in terms of a boson model. Such an ap- 
proach, called bosonization (see, e.g. [1,0|), is well known 
and successful for ID fermionic systems. Attempts to 
bosonize fermionic models in the dimensionality > 1 
have been undertaken in the past, starting from the works 
lol . 11 1 followed by [l^l- These schemes, however, have 
problems when dealing with large momentum transfer by 
the interaction [ist . They reproduce the random phase 
approximation (RPA) but do not reach beyond. 

A more general low energy bosonization scheme sug- 
gested recently [1^ is based on quasiclassical equations 
and takes into account large momentum transfer. New 
logarithmic contributions to anomalous dependence of 
the specific heat Q and spin susceptibility [l^l were 
found. However, working well for d = 1 the scheme of 
Ref. 1J| is not completely accurate for c? > 1 missing 



hence they cannot be used for accurate numerical studies 
of the initial fermionic problem. 

In this paper, we present a new scheme that allows 
one to map interacting fermions to interacting bosons 
exactly. This mapping works in any dimension at any 
temperature. The effective model obtained describes in- 
teracting bosonic excitations. It can be written either in a 
form of a model of non-interacting bosons in a Hubbard- 
Stratonovich (HS) field with a subsequent integration 
over this field or in a form of a field theory containing 
superfields with quartic and cubic interactions. The for- 
mer may be convenient for MC study, while the latter 
promises to be good for analytical investigations. 

We start with a general model of interacting electrons 
described by the Hamiltonian 



H — Ho + Hint, 

where Hq is the bare part, 

and Hinf stands for an electron-electron interaction. 
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In Eqs. (Pis]) Crcr (c^ct) ^^'^ annihilation (creation) opera- 
tors of the electrons on a lattice site r with spin cr = ±. 
The function t^-.r' describes the tunneling from the site 
r to the site r', T4,r' is the electron-electron interaction 
between the r and r' and fi is the chemical potential. 

The scheme of the bosonization suggested here can be 
developed for arbitrary functions trr' and Vrr' in an ar- 
bitrary dimension. However, in order to make formulas 
more compact we assume that 



some effects of the Fermi surface curvature [161 



All the previous bosonization methods are not exact; 



Vr,r' = SryVo, Vq > 0, 

which corresponds to an onsite repulsion. 



(4) 
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Then, the term Hi„t can be rewritten in the form 



while replacing the chemical potential /i by /i' = ^— Vo/2. 

In this paper, we concentrate on studying thermody- 
namics and calculate the partition function 



(6) 



Z = Ttcxp{-(3H] , (3^1/T. 



As Hq and H^^^l do not commute we subdivide the in- 
terval (0,/?) into slices of the length A = j3/N <C (3 and 
write Z as a time ordered product over the imaginary 
time r. Following a standard route of the HS transfor- 
mation we decouple the interaction term integrating 
over a periodic real field 0,. (t) = 0,, (r + /?) and come to 



Z = 



TV 



hm / ^ [0] exp -^Y.Y. ^r-i (7) 
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Cr^a,i and Cr,cr.i are the annihilation and creation oper- 
ators in the interaction representation[lj taken at t; = 
{I — 1/2) A, D4) is the normalized product of all differ- 
entials d(pr,i and (pr.i = ipr {ti). The product in ^ [(/)] is 
ordered in time, such that Z = 1 is on the right. 

For the analytical study, we could write Eq. ([7]) ex- 
plicitly in the continuous limit A — > using integrals 
and time ordering operators T^. However, MC calcula- 
tions imply finite A with typical values of (jy^j growing as 
A~^/^ for A ^ and therefore, we keep finite A. 

Calculation of the trace over the fcrmionic operators 
in Eq. ([7]) is not simple for finite A and one should 
approximate Z [0] by a more convenient expression. A 
standard approximation Zf [</>] used in MC simulations 
instead of Z [cj)] , Eq. reads [i,[6| 



N 



det 



l + ncxp(-/i[0(TO]A , (8) 
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K.a [0 (t)] = ir ~ - Cr0r (t) 

where irfr = — X^r' ^r.r' fr' for an arbitrary function fr- 
We suggest here another approximation Zb [0] to Z [(/)] 
that is more suitable for the bosonization. 



Zb[(f'] = Tr 



cxp 



Ho 
T 



(9) 



exp (| ^ I Cr(f>r (r) Cr^a (t) Cr,a (t) cIt 



0^ (r) = for - 1) A < r < ^A. 

The functional Zb[(j)], Eq. (O, differs from Z [(/)], Eq. 
([7]). by integration of the operator c^.o- (r) c^^cr (t) over 



each slice instead of taking it in the middle of the slice 
and multiplying by A. Therefore, the difference between 
Zb [(/)] and Z [(/)] should vanish in the limit A 0. We 
emphasize, however, that Z[(j)\, Zf [(/>], and Zb [4>\ differ 
from each other at finite A. The functional Zb [(/>] , Eq. 
([5]), has a form of the exact partition function for an 
electron in an external (generally, discontinuous in time) 
field (/),■ (t) and we can use standard transformations. 

In order to reduce the fermionic model, Eqs. ([9]), to a 
bosonic one, we introduce as in Ref. [l4| an additional 
variable < u < 1 and write the function Zb [4>\ as 



Zb[cl)]^Zocxp[Y.l^ I 



(J4>r (t) Gi'^^^l (t, T + 0) dudT 



where Zq is the partition function of the ideal Fermi gas 
and G^"^-*^ (''";''"') is a fermionic Green function, 



d 



with the boundary conditions 



M{t) ]G'^;,l{Ty)^6ry5{T-T'), 

(10) 



-G 



(r + /3,r') 



We develop our bosonization scheme introducing 



Ary {z) = d^l, (r, r + 0) - G^^', (r, r -f 0) 



(U0) 



(11) 



where z = {T,a,u) and g|.°^, (r, t') is the bare elec- 
tron Green function. The function Ar,r' (t) is periodic, 
Ary (t) = Ar^r' {t + P) , and, hence, describes bosons. 
We rewrite the partition function Z [(jj] as 



Zb [0] = Zq cxp - 



'- „ - Jo Jo 



4>r (t) Ar^r {z) dudr 



(12) 

and derive a closed equation for Ary (z)- For that pur- 
pose we write a conjugated equation 

A _ h,.,^^ ^u4> (r')] ) Gl!:Sl {r, r') = S^yS (r - r') 

and subtract it from Eq. (jlOp putting in the resulting 
equation r' = t + 0. The same can be done for the bare 
Green function g|."^,.^ (r, r') and finally we obtain 



d 

— — j4j. r' {z) + T~Lr.r' (,z) Ar^r' {z) 

cIt 



-ua^r,r' (t) riry , 

[z) = ir — £r' — Ua^r.r' (t) , (13) 



^ry (t) = 0r (t) - 0r' (t) 
.(0) 



The function n^.^' = G; (r, r + 0) in Eq. ^ is the 
Fourier transform in r ~ r' of the Fermi distribution. 
Equation ((T3|) should be supplemented by the condition 



^Ar^r{z) = 0. 



(14) 
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Equation ((T4|) can be obtained noticing that Ar^r (z) 
is a constant independent of (f>r (t) , which follows from 
Eq. (|13p . Assuming that the interaction and, hence, 

(r) vanishes at infinity we come to Eq. 

So, we are to solve the linear equation (fT^ for Ar^,.' (z) 
with the condition (HH), substitute the solution into Eq. 
(fT2)) and then into the first equation ([7]). A possible 
strong discontinuity of the function </> does not lead to 
any problems in the limit A ^ 0. Our scheme is sim- 
ilar to that of Ref. [3| developed in the quasiclassical 
approximation but now all the transformations are exact. 

It is convenient to exactly integrate over the field (f)r (r) 
in the beginning and thus derive a field theory for inter- 
acting bosons. In Ref. [IJ] this goal has been achieved by 
integrating over 48-component supervectors, which has 
led to a rather cumbersome Lagrangian. Now we use 
another trick, known as the Becchi-Rouet-Stora-Tuytin 
(BRST) transformation, based on introducing superfields 
[l3| (see also the book A similar transformation 

was used in the quantization of non-abelian gauge theo- 
ries [3| • In condensed matter physics this trick has been 
used first in Ref. 

Within this method one replaces solving an equation 



F{A)^0, 



(15) 



where is a matrix function of a matrix function A, and 
a subsequent calculation of a quantity B{Ao), where Aq 
is the solution of Eq. (|15p . by an integral of the form 



det ( — 

oa 



da. 



(16) 



B = J B{a) 6 [F {a 

The 5-function can be written as 

S[F{a)]=C J exp[ifF{a)]df, 

where C is a coefficient, and the determinant is obtained 
after integration of an exponential of a quadratic form in 
Grassmann variables rj and 77+. 

Our problem of solving Eq. and calculation of the 
integral in Eq. is of this type and we proceed follow- 
ing the above trick. We introduce anticommuting vari- 
ables 9 and 0* and a superfield *r,r' (R), R = {z,0,O*}, 

^r^r' (R) - ar^r' (z) 9 + f^,. (z) 9* + T^^y {z) + r,+y {z) 9*9 

where a, / are real and rj is an anticommuting field. The 
field 5* is periodic as a function of r, [t) = [t + (i) , 
but is anticommuting. The hermitian conjugation " -f " 
implies both the complex conjugation and transposition 
"T" with respect to r, r'. 

As a result, one comes to an effective action quadratic 
in 5* and linear in 4>r{T). This allows us to integrate over 
(t)r (r) with the Gaussian weight of Eq. ([T]) and we come 
to the final expression for the partition function Z , 



Z^ZoJ exp [So [*] - Sb - Si [^]] D^, 



(17) 



where Sq ['^] is the bare part of the action, 

d 



So 



2^ J r' \dT 

and the interaction terms are given by 



dR, 



^r.r{Rl)0l+2lIlr{Rl) 



aaidRdRi 



Si 



2 ^ 



S{T~Ti)Ilr (R) lir (i?i) aaidRdRi, 



Ur (i?) = It ^ [ (^-r'.r (R) - Tl,.' ,rO) (*r,r' (R) - Tir 



Integration over R in Eq. (jl7[) implies summation over a 
and integration over u,t,9,9*. The bare action So and 
the interaction term 5*/ are invariant under the transfor- 
mation of the fields 'J 



9,9*) '^ry (61 4- 



(18) 



with K and k* being anticommuting variables, whereas 
the term Sb breaks the invariance. The invariance un- 
der the transformation (fT8|) is stronger than the standard 
BRST symmetry for stochastic field equations (invari- 
ance under the transformation 4* (6**) {9* + k*)), 
Ref. , and reflects additional symmetries of Eq. (jl3p . 
It differs from the full supersymmetry by the presence of 
the term nrir^r' in Eq. psp but still can lead to interest- 
ing Ward identities. 

The model described by Eqs. pT)) can be studied using 
standard methods of field theory. One can, e.g., expand 
in the interaction Vq or develop a renormalization group 
scheme analogous to that of Ref. [l3| ■ In both the cases 
one can use the Wick theorem with rather simple contrac- 
tion rules that follow from the form of the bare action So . 
We leave such calculations for future publications. 

Neglecting cubic and quartic in 5* terms in Sb and Si 
in Eqs. p7p one has a purely quadratic action and the 
partition function Z yields an RPA-like expression. 



Z Zo exp 



K^l + Vo 



(19) 



2^^-^ (2^) 

7i(p-k/2) -n(p + k/2) d'^p 
zu; + e (p - k/2) - £ (p + k/2) (2^ f 



The same result can be obtained using Eqs. p^ [T3)) and 
neglecting the field ^^(t) in the L.H.S. of Eq. p3|l . 

In Eq. p9|) . (A' — 1) is the contribution of non- 
interacting bosonic excitations. Considering their inter- 
action originating from the cubic and quartic in term 
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in Eqs. (|17p one can fully describe the initial fermionic 
system. So, going beyond RPA, Eq. (fT9|) . is straightfor- 
ward and this is a very important advantage with respect 
to the older bosonization schemes [13, U, 12, We 
are confident that the present scheme can improve the 
analysis of Ref. of non-analytical corrections to the 
Landau Fermi liquid theory and expect its usefulness for 
study of a large variety of problems of strongly correlated 
systems. 

Now we sketch a possible route for MC simulations. 
Standard MC algorithms are based on using Eq. ([5]). 
However, for some important configurations of 0^ (t) the 
fermionic determinant Zf [(j) (r)] is negative, which makes 
the MC procedure inefficient. This is the famous sign 
problem. 

Here we suggest to use instead of Zj [0 (r)] the func- 
tional Zi,[(j){T)], Eq. p2)) . that can be found solving 
Eq. (HSl). The solution of Eq. ^ and the function 
Zb [(j) (r)] can be approximated using a Green function 
Qr'r''^ri r' (''"' '''i ) iutroduccd as the solution of equation 

Then, we write the functional Z^ [(j) (r)] as 



Zb [(t> (t)] = Zo exp 



- E E 

cr,r,ri,r[ i.j 



(20) 



Similarly to Eq. 



we write the function 



^r.'rTi,'-; '^i) > > Tj > in the form 

g:::,Uy^{n,r,)=P^y{n,r,) (21) 

X ^1 - Pry (Tj-,0) A-,r' (/3, Tj)) Sr^riSr',r{- 

Herein, the operator P is given by the expression 

Pry in, Tj) = Yl exp i-Hry (r,, a, u) A) (22) 

i>l>3 

where the multipliers in the product are ordered in time 
growing from the right to the left. (Of course, one 
should discretize also the variable u). The function 
g'^'^,'^^^ ^, {Ti,Tj) satisfies the symmetry relation 



(r, T-) = _Q'^;^<P 



(23) 



which allows one to consider times (3 > tj > Ti > 0. 

The form of the Green function Eq. ((2T|) is typical for 
bosons. By construction (see Eqs. (fT3l [2T |l ) it is real un- 
less a singularity is present, in which case an imaginary 
part might be generated. We argue that a possible zero 
in the Bosc-dcnominator in Eq. (pij) is compensated by 



the function y , Eq. (fT3|) , vanishing at ri = r[. Al- 
ternatively, one can antisymmetrize in the beginning the 
function Q in ri, r'l by antisymmetrizing the (5- functions 
in Eq. (|13p . This compensation is clearly seen in the 
RPA, Eq. HH). 

In the absence of any singularity, the result is insen- 
sitive to the way of subdividing the interval (0, f3) into 
slices and Zb[4> (t)] remains positive in the process of the 
calculation for any (j)r (t^) excluding the sign problem. 
Since Z can now be expanded in a sum of positive terms, 
we believe that this MC procedure can be efficient. The 
above derivation can be done using the "Ising spin" auxil- 
iary field of Refs. [1,0] as well, which is usually preferable 
for MC computations. 

In conclusion, the exact bosonization method pre- 
sented here opens new possibilities of both numerical 
and analytical study of models of interacting fermions. 
There is a reasonable chance that this new formalism is 
free from the sign problem supposed to be generically 
NP-hard or problems of equivalent complexity. 
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